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CONSTRAINTS ON SHALLOW 56 NI FROM THE EARLY LIGHTCURVES OF TYPE IA SUPERNOVAE 
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ABSTRACT 

Ongoing transient surveys are presenting an unprecedented account of the rising lightcurves of Type la 
supernovae (SNe la). This early emission probes the shallowest layers of the exploding white dwarf, which can 
provide constraints on the progenitor star and the properties of the explosive burning. We use semi-analytic 
models of radioactively-powered rising lightcurves to analyze these observations. As we have summarized 
in previous work, the main limiting factor in determining the surface distribution of 56 Ni is the lack of an 
unambiguously identified time of explosion, as would be provided by detection of shock breakout or shock- 
heated cooling. Without this the SN may in principle exhibit a "dark phase" for a few hours to days, where 
the only emission is from shock-heated cooling that is too dim to be detected. Nevertheless, by considering 
the time-dependent velocity evolution, the explosion time can be better constrained, albeit with considerable 
uncertainty. This technique is used to infer the surface 56 Ni distribution of three recent SNe la that were caught 
especially early in their rise. Although we cannot constrain the explosion times to better than « 1 day, in all 
three we find fairly similar 56 Ni distributions. Observations of SN 201 lfe and SN 2012cg probe shallower 
depths than in SN 2009ig, and in these two cases 56 Ni is present merely w 1O~ 2 M from the WD's surface. 
We also use our conclusions about the explosion times to reassess radius constraints for the progenitor of 
SN 201 lfe, as well as discuss the roughly f 2 power law that is inferred for many observed rising lightcurves. 
Subject headings: hydrodynamics — shock waves — supernovae: general — white dwarfs 



1. INTRODUCTION 

Type la supernovae (SNe la) play a central role in modern 
astrophysics. They are use d as distance indicators to probe the 
expan sion of the Universe dRiess et al]l!998t iPerlmutter et alJ 
119991) . t hey produce most of the iron-group elements in the 
cosmos dlwamoto et alJ 1999b . and t hey provide an astrophys- 
ical c ontext for studying explosions dHillebrandt & Niemeverl 
120001) . But their importance has brought attention to the the- 
oretical uncertainties that frustratingly remain. It is gener- 
ally accepted that they result from unstable thermon uclear 
ignition of degenerate matter dHovle & Fowled [i960!) in a 
C/O white dwarf (WD), but the progenitor systems have not 
been identified. Candidates inclu de stable accretion fro m a 
non-degenerate binary comp anion dWhelan & Iberi 19731). the 
merg ing of two C/O WDs (llben & Tutukovlll984t iWebbinkl 
119841) . or accretion and detonation of a helium shell on a 
C/O WD that leads to core detonation dWooslev & Weaver] 
Il994at ILlvne & Arnettlll995l) . In addition, it is not known 
whet her the incineration procee ds as a sub-sonic deflagra- 
tion dNomoto et al.l 119761 119841) or deflagration-detonatio n 
3vll99ll: F" 



becoming more common, the time is ripe to explore what can 
be learned from these measurements. 

In the following work we use semi-analytic models to 
study where and how much 56 Ni is present in the outer 
ejecta of SNe la. As discussed in ou r previous investigat ion 
of radioactively-powered lightcurves dPiro & Nak ar 2012), it 
is difficult to directly measure the 56 Ni distribution with- 
out a detection of the explosion time, as woul d be pro- 
vided by shock breakout or shock-heated cooling dPiro et al 



2010 tlNakar & SarilHoToHNakar & Sarill^HllRabinaketal 
20121) . Unfortunately, in the case of SNe la, such emission 



transition (DDT: lKhokhlovll99UlWooslev & Weaverll994bl) . 
Each scenario has implications for the velocity profile, density 
structure, and distribution of ashes within the exploding WD. 

A powerful method for constraining between these mod- 
els is the study of the early-time behavior of SNe la, since 
this is when the shallowest layers of the WD are probed by 
the observed emission. Analysis of spectra provides one way 
of le arning about the surfac e abundances of these explosions 
(e.g.. lHachinger et al.l 120121) . The photometry is also s ensi- 
tive to the depth and distr ibution of radioactive heating dPirol 
l2012tlPiro & Nakaril2012l) . With early observations of SNe la 
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has never been detected because of the small WD radius. If 
merely photometric lightcurves of the rise are available, there 
is a degeneracy between emission being from 56 Ni near the 
surface with a recent explosion versus 56 Ni deeper in the star 
but with an explosion further in the past. In the latter case, a 
SN la exhibits a "dark phase" for a few hours to days until the 
thermal diffusion wave reaches the shallowest 56 Ni deposits. 
Even with these uncertainties, constraints can still be provided 
by comparing a wider range of properties, such as the velocity 
evolution. This information is available for a few well-studied 
SNe la, and we use it in order to estimate the time of explosion 
and surface 56 Ni distribution for each of them. 

In §|2] we summarize the semi-analytic framework used to 
model the rising lightcurves. In §|3] we analyze observations 
of three recent SNe la and summarize our constraints on their 
shallow 56 Ni distributions. In ^4]we consider the t 2 rise that is 
often observed in early lightcurves and discuss whether t 2 (or 
any power-law) should be expected. We conclude in ^5] with 
a summary of our results and a discussion of future work. 

2. RADIOACTIVELY-POWERED RISING 
LIGHTCURVES 

In the following we present the model used for this study, 
which borrows from and builds upon our recent work on 
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radioactively-powered rising lightcurves. In iPirol (120121) we 
focused on direct 56 Ni heating at the depth of the diffusion 
wave. In lPiro&Naka^dloTS we added the "diffusive tail," 
which provides heating at depths shallower than the intrinsic 
56 Ni distribution. Here we include these effects in greater de- 
tail by integrating over their co ntributions throughou t the WD, 
as discussed in Appendix B of lPiro &NakariJ20T2h . 

As the ejecta from the SN expands, a thermal diffusion 
wave travels back through the material. At any time t, the 
diffusion wave has a depth of roughly 
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where E = 10 51 £'5ierg is the explosion energy, and M = 
\AM\aMq is the ejecta mass, and k = O.Iko.i cm 2 g" 1 
is the opacity. We approximate the opacity as constant, 
which is roughly corre ct for the temperatures of interest 
dPinto & Eastmanl2000l) . T he scalings and prefactors in equa- 
tion CQ) use Appendix C of iPiro & Nakari (120121) with values 
appropriate for Chandrasekhar mass WDs. 

At early times the ejecta is optically thick to gamma-rays 
emitted from radioactive decay, and they efficiently heat the 
SN. Heating in material shallower than AMdiff directly goes 
into the observed luminosity. Heating in material deeper than 
AMdiff only contributes to the observed lightcurve if some 
fraction of the photon s from these larger de pths are able to dif- 
fuse up to AMdiff. In lPiro & Nakari d2012l) . we demonstrated 
that if some depth in the ejecta would be reached by the diffu- 
sion wave at time t', then at times t <t' the fraction of photons 
that would escape from this depth are 

Escaping fraction w erfc(r'/\/2f) (2) 

where erfc is the complementary error function. This pro- 
duces the so-called diffusive tail. 

Motivated by this picture, we split the total observed lumi- 
nosity into two parts 

L(t) = L 56 (t)+L tail (t), (3) 

where L56 is the direct heating by 56 Ni down to AMdiff, and 
Ltaii is the diffusive tail from material deeper than AMdiff. 
Each is an integral over different regions of the ejecta. For 
the direct heating component 

L 56 (t) = j'x 56 (t')^^e(t)dt', (4) 

where Xs(,(t) is the mass fraction of 56 Ni at the depth of the 
diffusion wave at time t , and the specific heating rate is 



e(0 = e Nl e-'/ fN - + eco(<T f/fc ° - e"' /fN '), 



(5) 



where e Ni = 3.9 x 10 10 erg g" 1 s" 1 , f Ni = 8.76days, e Co = 7.0 x 
10 9 erg g" 1 s -1 , and fco = 11 1.5 days. The total diffusive tail 
component is the integral over all the diffusive tails from heat- 
ing deeper than AMdiff, 



L t ail(0 : 
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dt' 



erfc(l/v / 2) 



We take the upper integration limit to be the diffusion time 
through the entire ejecta fdifO. which roughly corresponds to 

3 Note that in this w ork we are using a different definition of fjjff than that 
in |Piro~& Nakar (20TJ). 



the time of lightcurve peak. 

Since AMdiff oc f 1 ' 76 , equations (0]i and © are rewritten as 



L 56 (f)=1.76Lo(f) 
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where Lo(0 = X5(,{t)AM^n{t)e{t) is roughly the local heating 
rate from 56 Ni. The luminosity has no contribution from the 
diffusive tail once the diffusion wave has travelled through the 
ejecta, thus we define Ldiff = Ls6(f = fdiff)- 

When actually performing calculations, it is useful to write 
these expressions in dimensionless forms. First, let x = t/tdis 
and x' = f'/ f diff, where x and x' vary from to 1. We define 
the ratio of the local heating rate to Ldiff as 

Mix) 



A(x)=1.76 
= eW 



Ldiff 



X 5b {x') fx 
X 56 (x) 
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The ratio of the observed time-dependent luminosity to Ldiff 
is then 



m 

Ldiff 
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In this form the right-hand side is dimensionless and only de- 
pends on the 56 Ni distribution. This allows us to vary Xs(,(x) 
and calculate a wide range of lightcurves, which can then be 
rescaled to a particular observation via Ldiff and fdiff- When 
fitting a 56 Ni distribution to a given lightcurve in the next sec- 
tion, we use the parametrization 



X 56 (x) = 



x; 6 



1 +exp \-(3{x-xx/2) 



(11) 



where X' 5b sets the normalization, /3 controls the steepness of 
the rise, and xi/ 2 is the time when Xse/X' i(t = 1/2. This al- 
lows us to consider a variety of 56 Ni distributions with two 
parameters. The normalization is determined by 

Ldiff 



4r 



1.76AM d iff(fdiff)e(fdiff) 
x' 0J6 dx' 



(i 



1 + exp [-/3(x' - 
and thus is not a free parameter. 



Xl/2) 
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3. COMPARISONS TO SPECIFIC SUPERNOVAE 

Recent observations have been especially fruitful in catch- 
ing SNe la at early times. We use this work to analyze three 
well-studied events: SN 201 lfe, SN 2012cg, and SN 2009ig. 
For each we summarize what can be constrained from their 
photometric lightcurves and velocity evolution. Although 
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there are particular issues for each event (which we discuss 
below), our general strategy is as follows. 

1 . Since a SN may in principle exhibit a dark phase, we 
assume that the time of explosion is not known. 

2. For a spectral line generated at constant specific opac- 
ity, it s velocity is a pow er law with time with v oc 
r° 22 dPiro &Nakarll2012l) . We vary the explosion time 
and check when the observed absorption features best 
match this power-lawQ. From this we infer what is the 
likely explosion time. 

3. The photospheric velocity v p h is expected to roughly 
follow the low-velo city Si II A6355 absorption feature 
dTanaka et al.l 120081) . Using the fits performed in the 
previous step, we can therefore estimate v p h(f)- The 
photospheric radius is then given by r p h = v p hf . 

4. Assuming that the SN emits roughly as a blackbody 
and using the observed B, V, and R lightcurves, we fit 
the color temperature T c and bolometric luminosity as a 
function of time using L w Airr^a^T^ . 

5. Theoretical lightcurves are generated with different 
X^{x) via equation ( [Tol l, where X^(x) has the func- 
tional form of equation ( TTTb . We estimate Ldiff and 
fdiff as roughly the peak luminosity and time of peak 
luminosity, respectively. In this way the theoretical 
lightcurves are rescaled for comparison with the bolo- 
metric lightcurve, and we can put constraints on what 
is the most likely distribution of 56 Ni. 

3.1. Modeling SN 201 lfe 

We first focus on SN 201 lfe because it is the most con- 
strained by our modeling. SN 201 lfe explod ed in August 
201 1 as the closest SNe la in the last 25 years dNugent et alj 
1201 lb . The considerable interest in this event and its prox- 
imity make it one of the best studied SNe la. The time- 
de pendent velocities o f absorption features are summarized 
in iParrent et alJ (120121). The B, V , and R rising lightcurves 
are presented in lVinko et al.l (120121) . This particular work was 
chosen because of the high density of observations during the 
rise, but we could have just as well considered other data sets 
dMunari et al .1120 12t iRichmond & Smithll2012l). We use a dis - 
tance modulus for M101 of 29.05 (IShappee & Stanekll201 ll) . 
and no reddenin g is included becau se it has been inferred to be 
relatively small (fPatat et al.l[201 ll) . The various studies have 
found different time for the peak bolometric luminosity, de- 
pending on the fitting method used. For the present work 
we choose a time of peak of JD 2455815.4, although this 
choice does not greatly impact our conclusions for the 56 Ni 
distribution at s hallow depths . The earliest detection was at 
JD 2455797.65 dNugent et al.ll2011l) . and their fitting of a t 2 
power law to the rising luminosity, (as is common practice) 
gives an explosion time of JD 2455797.2 before the peak. 
Although they quote an error of ±0.01 days, as we discuss 

4 The fact that we find that all of the absorption features obey roughly 
the same power-law dependence is interesting in and of itself and deserves 
some discussion. It indicates that different features are due to different line 
opacities within a flow with the same velocity power-law profile and are not 
separate velocity components in the ejecta. This argues against situations 
where the high velocity fe atures are g enerated by a separate event during or 
prior to the explosion (e.g,. IPirrJ201 11) . 
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FIG. 1 . — The top panel plots the \ 2 found by fitting the velocity evolution 
as a power law with time for different explosion times and various power- 
law indices as labeled. We draw vertical dot-dashed line at the time of an 
upper limit from a non-detection and a shaded region at the explosion time 
inferred by fitting a r 2 rise ( Nugent et al. 201 1). In the bottom panel we plot 
the observed low-velocity Si II A6355 (circles), high-velocity Si II A6355 
(squares), and high-velocity Ca II H&K (crosses). Filled and open symbols 
indicate data that was used or not used for the fit, respectively. The lines show 
our best fit velocity evolution for v oc f ~° 22 , and the solid line indicates the 
Vph we use in subsequent analysis. 

later this practice is not well justified and the tr ue uncertainty 
in the explosion time is considerably larger. iNugent et alJ 
(1201 ll) also observed the location of SN 201 lfe roughly at 
JD 2455796.7, which provides an upper limit in the apparent 
g-band magnitude of 21.5 (absolute magnitude of -7.55). 

In our analysis of the velocity evolution, we use the low- 
velocity Si II A6355, high-velocity Si II A6355, and high- 
velocity Ca II H&K absorption featurefl In the top panel 
of Figure Q] we plot the x 2 found by fitting these features 
with power-law velocity profiles as a function of the explo- 
sion time, where \ 2 is defined as 

2 sr^ f v N -v(t)\ 2 
N \ / 

where N is the number of data points, vjv is a measured ve- 
locity, and Av = 500 km s -1 is a rough estimate of the mea- 
surement error (Par rent et al.ll2012T) . The reduced x 2 around 
the best fit explosion time is about 1.2 (there are 17 degrees 
of freedom). In Figure Q]we consider three different velocity 
power-law indices centered around the model prediction of 
v oc f -0 - 22 . This shows that the model provides a good descrip- 
tion of the data, and that assuming that the power-law index is 
known, the explosion time is measured to within about ±0.25 

5 Other absorption features are measured, but we restrict our study to these 
three since they are some of the most widely available in SN la literature. 
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FIG. 2. — A summary of the fits to SN 201 lfe. The top panel shows the 
inferred bolometric luminosity (filled circles) and fit bolometric luminos- 
ity (solid curve), which is composed of L56 (dashed curve) and Ltail (dotted 
curve). The middle panel shows the color temperature, and the bottom panel 
shows M56 = Lss/e. 

days. However, assuming slightly different power-laws pro- 
duces fits with similar quality and results in explosion times 
that vary by about ps 1 day. Since theoretically v oc t~° 22 is 
the preferred velocity profile we consider JD 2455796.6 to be 
the most likely explosion time with an uncertainty of roughly 
±0.5 day. This is actually very simila r (within 0. 1 days) of the 
non-detection bv lNugent et al.l d201 ll) . In the bottom panel we 
present the velocity data along with our best-fit velocity evo- 
lutions. Open symbols indicate data that were not used for the 
fit because they are near peak where the velocity profile is not 
expected to be a power law. 

For any given explosion time we can look for the 56 Ni that 
produces the observed luminosity. The results from fitting the 
photometric observations are presented in Figure [2] In this 
particular case we use the time of the non-detection for the ex- 
plosion time, which is sufficiently close to our preferred time 
that the qualitative features are unchanged. In the top panel 
we compare the inferred bolometric lightcurve (filled circles) 
to the model fit (solid curve). We also plot the contributions 
from direct heating L56 (dashed curve) and the diffusive tail 
Ltau (dotted curve). The direct heating component is larger at 
late times and the diffusive tail is stronger at early times when 
the 56 Ni is less abundant. Nevertheless, L ta ji is never more 
than a factor of two greater than Lg6. This means that at least 
very roughly, the observed bolometric luminosity can be used 
to infer the distribution of 56 Ni, and that 56 Ni must be present, 
at least in some amount, at the depths that are probed by the 
earliest emission. 

The middle panel of Figure|2]shows the inferred color tem- 
perature T c , and the bottom panel shows the mass of 56 Ni 
above the diffusion wave depth, given by 



M 56 (t) = L 56 (t)/e(t). 



(14) 
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This is roughly independent of the explosion time because it 
is just set by the bolometric luminosity at any given time. In 



Time after explosion (days) 

FIG. 3. — Comparison o f the early g-band data I Nugent et al. 2011]) and a 
non-detection upper limit (Bloom et al. 2012) to theoretical lightcurves from 
radioactive heating (da shed curves) and shock-heated cooling (solid curves) 
calculated according to Piro et al. (2010). This does not include the suppres- 
sion of the shock-heated cooling (or "drop out") that occurs when the dif - 
fusion wave moves into ideal gas dominated material (Rabinak et 311 120121) . 
The top panel is roughly the explosion time inferred from a t 2 extrapolation. 
The bottom panel assumes that the explosion occurred 0.5 days earlier, for 
which the radius constraint is a factor of 1.9 larger. 

contrast, T c changes with explosion time because an explo- 
sion further in the past implies more expansion at any given 
time and thus a smaller T c . This means that an additional con- 
straint on the explosion time could be made via a temperature 
measurement, although this requires detailed spectral model- 
ing that i s outside the scope of this work (see the discussion 
of ^ in lPiro & Nak^l2012h . 

3.2. Radius Constraints and Shallowest 56 Nifor SN 201 lfe 

Using the data from iNugent et al.l d201 ll) and a non- 
detection rj 7hrs earlier. IBloom et all (120121) argued that 
the progenitor of SN 201 lfe had a radi us < 0.02K W by 
using models of s hock-heated cooling dPiro et al.l 120101 : 
iRabinak et al.ll2012l) . But this assumed that the time of ex- 
plosion could be accurately dete rmined from extrapol ating t 2 
back in time. As emphasized in[Piro & Nakai] d2012l) . this is 
not generally a robust method for finding the explosion time 
(see also @, so it is worth revisiting the radius constraint for 
a range of explosion times. 

In Figure [3] we plot the early data and non-detection upper 
limit for SN 201 lfe for two different explosion times. The 
theoretical curves include radioactive heating (dashed curves) 
and shock-heated cooling (solid curves). The first thing to 
note is that 56 Ni cannot always be present at the earliest times 
and still produce the observed lightcurves. In the bottom 
panel we had to cut off the 56 Ni for times earlier than 0.9 days 
after explosion i n order to not overp redict the g-band upper 
limit reported in IBloom et al.l d2012l) . (In the top panel no 
56 Ni cut-off is needed.) This implies that for earlier explosion 
times there is a sharp cut-off in the 56 Ni distribution near the 
depth that generates the luminosity of the first detected light. 
This is not unexpected since 56 Ni probably does not extend to 
the very surface a nd th e earliest emission will be due to the 
diffusive tail. In 33 .51 we further discuss what depth in the 
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FlG. 4. — The same as Figure □ but for SN 2012cg. The shaded region 
shows the inferred explosion time from Silverman et al. 12012) using r. 

WD is implied by this time. 

The other thing to note from Figure [3] is that when the ex- 
plosion time is further in the past, upper limits on the emission 
from shock-heated cooli ng (solid curves) are not as stringent. 
Using the models from iPiro et all ( 120101) we find that when 
the explosion is merely 0.5 days further in the past (the bot- 
tom panel) the radius can be a factor of 1 .9 greater than in the 
top panel. 

Another potentially important effect that is not included in 
Figure [3] is the "drop out" in the shock-heated cooling emis- 
sion that is expected once the diffusion wave exposes the 
depth where the shock is matter rather than radiation domi- 
nated. This is expect ed to occur ~ hours after explosion for 
a typical WD radius dRabinak et alJl20T2h . Although we do 
no t consider explosion times earlier than the non-detection 
of iNugent et al.l (1201 ll) in Figure [5] it is possible that the ex- 
plosion occurred before (~ 0.5 day ) this time, and the non- 
detection is simply during the dark phase between the drop 
out and the latter 56 Ni heating. iBloom et ail d2012l) find that 
the drop out limits the radius constraint posed by their upper 
limit. If the explosio n is one day before the date estimated by 
INugent et al.l (1201 ll) . then the uncertainty in the limit on the 
radius of the progenitor of SN 2010fe is only ~ 0. IRq- 

3.3. SN2012cg 

The velo cities and photometry for SN 2012cg are sum- 
marized in [Silverman et al j20 12). Further photometry is 
presented by iMunari et al.l ( 20121) . including data around the 
peak identified to occur at roughly JD 2456083.0. The ve- 
locity fitting results are shown Figure [4] (again taking Av = 
500 km s" 1 ). The best fit explosion time is JD 2456063.5, 
but the strength of the fit is not as strong as for SN 
201 lfe. In comparison, using t 2 ISilverman et al.l ( 120121) find 
JD 2456063.2 ±0.2 (the shaded region in the top panel of 
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FIG. 5. — The same as Figureff] but for SN 2012cg. 

Figure [4), which is consistent with our fits. The lightcurve 
modeling from the observed B, V, and R measurements use a 
distance modulus of 30.9. The summary of our results from 
the photometric data are presented in Figure [5] 

3.4. SN2009ig 

The v elocities and photo metry for SN 2009ig are pre- 
sented in iFolev et alj (120121) . The time of B-band peak is at 
JD 2455080.54, and the distance modulus is 32.6. The evolu- 
tion of the Si II A6355 absorption feature is a little more com- 
plicated in this case and deserves some discussion. At early 
times (earlier than 12 days before fi-band peak), Si II appears 
to only have a high velocity component, and a low velocity 
component grows to be more prominent later. We take the 
low velocity component as indicative of the photosphere, but 
use both the high and low velocity components when fitting 
the v oc f "° 22 power law. Data taken when the features overlap 
could potentially bias the fit due to blending, but we did not 
find that it has an adverse impact on our fits. 

In Figure [6] we summarize the velocity fitting. Only high 
and low velocity Si II are used in this case. High-velocity 
Ca II H&K absorption features may be blended wit h Si II 
A4130, and are not presented by IFolev et al.l (120121) . The 
best fit time of explosion is at JD 2455061.8. In compari- 
son, using a t 2 rise IFolev et alj d2012l) infer an explosion time 
JD 2455063.4 ± 0.07. Although SN 2009ig has the least con- 
straining fits of any of the SNe, this later explosion time seems 
difficult to reconcile with the velocity evolution unless v(t ) is a 
much shallower power law with time than that expected from 
theory. In Figure|7]we plot the best fit lightcurve properties as 
was done for the other SNe. 

3.5. Comparing and Contrasting Events 

In Figure[8]we plot the distributions of 56 Ni inferred for the 
three SNe la modeled above. In each case multiple values 
for the explosion time are considered to demonstrate how in- 
ferences on X56 change with this parameter. The solid lines 
in each panel indicate the preferred explosion time. For SN 
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FIG. 7.— The same as FigurefJ but for SN 2009ig. 

201 lfe (top panel), thick lines show the distribut i on cov ered 
by the photometric observations of IVinko et al.1 d2012f) and 
thin lines sho w the distribution in ferred by the earlier ob- 
servations by Nugent et al.l ( 1201 lb . This shows that since 
AMdiff oc f 176 , having observations just a day or two ear- 
lier can probe much shallower regions of the ejecta. For 
the preferred explosion time, X56 ~3x 10~ 2 at a depth of 



AMdiff ~ 10~ 2 M Q . T hese r esults are roughly consistent with 
models presented by iPirol ( 120121) . which assumed a similar 
explosion time but did not include the diffusive tail. As dis- 
cussed in 33.21 the upper limit on the luminosity at early times 
implies that there must be a cut-off in the 56 Ni distribution for 
some explosion times. These shallowest 56 Ni depths are indi- 
cated by filled circles in the top panel of Figure [8] (although 
not mentioned in 33.21 for the -0.5 days curve, 56 Ni cannot be 
shallower than 1 .7 days after the explosion). 

The 56 Ni distributions in SN 2012cg and SN 2009ig are 
fairly similar to SN 201 lfe over similar depths. The main 
difference is that SN 2012cg shows somewhat more 56 Ni 
around a range of AMdiff ~ 10~ 2 - 10 _1 Mq. Does this imply 
that SN 2012cg has more shallow burning products? Analy- 
sis of the spectra indicate that SN 201 lfe has considerably 
more unburned car bon at shallow depths than SN 2012cg 
dParren t et at 2012), which is at least consistent with this hy- 
pothesis. 

SN 2009ig also has a number of differences that are worth 
discussing. The Si II velocities at w lOdays past explosion are 
considerably higher in this event than either SN 201 lfe or SN 
2012cg. If this indicates a difference in the actual explosion 
energy, then using v p h oc E 039 dPiro &Nakaril2012l) argues 
that SN 2009ig was a factor of « 2 more energetic than the 
other two events. Such an explanation seems difficult to rec- 
oncile with the peak luminosity of SN 2009ig, which is fairly 
standard for SNe la. Another attractive possibility is that the 
large velocities are due to an asym metric explosion tha t is di- 
rected more toward the observer (Maedaet at 2010a|). For 
such larger velocities, there is more expansion and a generally 
cooler SN, as can be seen by t he T r presented in the middle 
panel of Figure □ iFoley et al.l (120121) note that SN 2009ig is 
considerably redder in the UV at early times in comparison 
to other SNe la and typical templates. Is this just due to the 
larger velocities? Another possibility is that these colors are 
due iron-peak elements near the surface, which again would 
be consistent with an explosion directed toward the observer. 
The mass fraction of 56 Ni for SN 2009ig is fairly similar to 
the other SNe at a depth of « 0. 1M Q , and data is not available 
early enough to probe shallower regions. 

Although the many differences found for SN 2009ig are 
tantalizing, we emphasize that these conclusions all hinge on 
our assumption that roughly v oc t~° 22 . If for some reason the 
velocity profile of SN 2009ig is different than the other two 
SNe, then these conclusions must be revised. On the other 
hand, if the velocity profile is significantly different in this 
case, that might be interesting in and of itself. 

3.6. Progenitor Models 

For all three SNe we study, 56 Ni must be present at least 
« 0.1M Q from the WD surface, and as shallow as « 1O~ 2 M0 
from the surface for SN 201 lfe and SN 2012cg (see Figure[8Jl. 
It is therefore worth discussing the implications for progenitor 
models and the chara cter of the explos i ve bur ning. 

As a comparison, lHachinger et al.l (120121) performed de- 
tailed UV/optical spectral modeling of SN 2010jn. From this 
analysis they also infer iron-group elements near the surface. 
DDT models can pr oduce 56 Ni near the WD surface (e.g., 
Ilwamoto et al.|[l999l) . but to get radioactive material as shal- 
low as « 10 ~ 2 M(7i may require a strongly mixed, off-center 
deflagration dMaeda et al .11201 Obi) . In DDT models with many 
ignition points that have fairly stratified ashes, radioactive ele- 
ments are not present near the surface. A gravitationally con- 
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FIG. 8. — Inferred distribution of 56 Ni as a function of depth in the WD. 
In each case we compare multiple explosion times, with the solid lines in- 
dicating the value preferred by fitting v oc r° 21 . The depth into the star is 
assumed to scale as AMdiff oc r 176 with a normalization of AMjiff = 1.4Mq 
at lightcurve peak. For SN 201 lfe (top pa nel), the thick curv es correspond 
to the constraints from the observations by Vinko et al. 12012), and the thin 
curves correspond to the observations by Nugent et al. 1 201 1). The filled cir- 
cles indicate the shallowest allowed deposits of 56 Ni so as not to overshoot 
the upper limit presented by Bloom et al. (2012). 

fined detonation also produces iron-pe ak elements near th e 
surface when a bubble rises and breaks dMeakin et al.ll2009l) . 

Another interesting scenario that may produce shallow ra- 
dioactive heating is the explosive ignition of a helium shell 
in a double-detonation. The depth and amount of 56 Ni we 
infer is not dissimilar to the helium shell masses needed for 
detonation and the total amount of radioactive material found 
for such events dShen & Bildstenll2009l:lFink et al.l20lol) . The 
main problem with such models is that if iron-peak elements 
are too abundant, they tend to produce colors that are too 
red and spectra that are inconsiste nt with normal SNe la 
dKromer et al.ll2010tlslm et al.ll2012l ). But if the helium burns 
in a lateral detonation which d oes not process the fu el as com- 
pletely to iron-peak elements dTownslev et al.l20"l2l) . this may 
overcome some of the difficulties double-detonation models 
have in reproducing observed SNe la. 

4. IS A t 2 RISE SPECIAL? 

A common practice with recent SN la observations is to de- 
termine the time of explosion by fitting the rising luminosity 
(often in a single band ) with a t 2 curve dNugent et ail 1201 It 
Milne & Brownl 120 12t iFolev et all 120121: [Silverman et al.l 



20121) . Studies of composite lightcurves formed from stack 



ing many SNe, which allow the power-law index to vary, find 
power-law indices of 1.8 ±0.2 dConlev et al.ll2006b . 1.8^;J| 
dHavden et al.l 120101) , and 2.20^ dGaneshalingam et al.l 
1201 11) . This then begs the question, is t 2 (or any other power- 
law) fundamental, and if not, what is the origin of these re- 
sults? 

Our discussion in ^2] shows that a priori a power-law lumi- 
nosity rise is not generally expected. The luminosity is driven 



by a combination of two factors: (i) the diffusion wave prop- 
agation, AMdiff(f), an d (ii) the distribution of 56 Ni fraction, 
X^. The expos ed mass does indeed evolve as a power law, 
with dPiroll2012l) 



AM diff (f ) oc t 



2(l+l/n)/(l+l/n+|8) 



(15) 



where n is the polytropic index and (3 is the power-la w index 
of the velocity gradient. For n = 3 and (3 = 0.186 (ISakurail 
119601) . this results in AM diff oc f 1 76 (as in eq. (TJ). In con- 
trast, the 56 Ni distribution is not well constrained by theory 
and may, in principle, vary in many ways. A power-law rise of 
the bolometric luminosity is expected only if the 56 Ni fraction 
evolves as a power-law as well, namely X$(, oc t a . In this case 
the bolometric luminosity evolves as L oc t L76+Q and since the 
photospheric radius is roug hly oc r 78 dPiro & Nakarll2012l) . 
the observed temperature ev olves rough ly as T c oc £(°- 2+Q 0/ 4 , 
This result was obtained by iPirol (120121) when the diffusive 
tail was not included, and we find that it still holds with the 
more detailed analysis presented in Since we do not ex- 
pect oc t °' 24 , our conclusion is that a t 2 rise (bolometric or 
in a single band) is probably not a generic property of SNe la. 
We also do not expect the rise to follow exactly any other 
power law. Moreover, since most explosion models predict a 
sharp decrease of in the outermost layers of the ejecta, the 
lightcurve is expected to rise exponentially (due to diffusive 
tail contribution) at very early times. How early this exponen- 
tial phase take place depends on the depth of the shallowest 
56 Ni deposit. 

What is then the explanation of the fact that analysis of large 
SNe samples are found to be consistent with a power-law rise 
with indices in the range « 1.8-2.2? It is probably a com- 
bination of two things. First, the unknown explosion time 
enables a reasonable fit even if the lightcurve is not exactly 
a power law. Second, in the depth range explored by most 
of these SNe rising phases, the is not varying by a large 
amount. This is because the first observation of most SNe la 
take place only a few days to a week after the explosion, so 
that the 56 Ni is distributed roughly uniformly or slowly in- 
creasing with depth. 

In Figure [9] we plot the bolometric lightcurves and fits 
for the three SNe we have been studying on a logarithmic 
scale to emphasize power-law dependencies. This shows that 
the lightcurves are not rising exactly as power laws, but that 
power law fits can provide a reasonable description of the data 
(although for SNe 200 lfe the rise is found to be slightly faster 
than f 2 ). This is because in all three of these SNe the X56 
is rising rather gradually over the depth range probed by the 
observations. 

There have been some attempts to explain why there should 
be a t 2 rise, but none of these provide arguments that are 
expected to hold in detail. The most simplified explanation 
is a fixed color temperature with a radius that increases lin- 
early with time ( Riess et al.ll 1 9991) . This model does not ex- 
plain why the temperature should be constant, and more im- 
portantly, in real SNe the color temperature does typically 
vary with time. A more fundamental explanatio n for a t 2 rise 
is give n in the Supplementary Information of iNugent et al.l 
d2011l) . which c onsiders radi oactive heating with thermal dif- 
fusion (also see lArnetHll982l) . This model makes two explicit 
approximations: (i) it ignores the velocity gradient, obtaining 
AM diff oc t 2 (basically setting j3 = in eq. HU), and (ii) it 
assumes that X^ is constant. Together these factors result in 
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FIG. 9. — The bolometric lightcurves for each of the three SNe (from the 
top panels of Figures [2] \5\ and [7), but in this case plotted with logarithmic 
axes to emphasize power-law behavior. 



a t 2 rise, but only for assumptions that are not realistic. 

To conclude, we expect the early rise to depend on the 
particular physical conditions in any given event and thus to 
possibly vary from one SN to another. It will be important 
to test this hypothesis in the future by building bolometric 
lightcurves from observations to infer just how much diver- 
sity there really is. Detailed numerical calculations of the 
rise will also be useful for understanding how much the early 
luminosity can change depending on composition and radia- 
tive transfer effects. Whatever the results are, extrapolating 
a lightcurve back in time with t 2 is not a reliable method for 
inferring the explosion time. 

5. CONCLUSIONS AND DISCUSSION 

Using early observations of three SNe la, and assuming that 
the absorption feature velocities evolve as v oc f" 022 , we con- 
strained the explosion times and shallow distributions of 56 Ni. 
We then used these findings to revisit the radius constraints on 
the progenitor of SN 201 lfe (in ^3. 21 , and discuss the t 2 rise 
that is reported for many SNe la (in $4} . Our general conclu- 
sion is that SN 201 lfe and SN 2012cg are very similar in most 



respects, including the rise time, 56 Ni distribution, and ener- 
getics. The main difference is that SN 2012cg has a slightly 
larger amount of shallow 56 Ni. 

SN 2009ig is somewhat different than the other two SNe. 
Although its 56 Ni distribution over the same depths probed in 
SN 201 lfe and SN 2012cg are fairly similar, it has higher ve- 
locities at any given time and its best fit time of explosion has 
the largest discrepancy with previous estimates (w 1.6 days 
earlier). This is curious because the peak luminosity of SN 
2009ig is fairly normal in comparison to the other SNe, and 
thus the amount of 56 Ni and the energetics should be simi- 
lar. One possible solution is if SN 2009ig is asy mmetric with 
higher velocities directed toward the observer (Maedaet al. 
I2010al) . Unfortunately our results on SN 2009ig are some- 
what tentative because it has the least constrained time of ex- 
plosion. This is because the low velocity Si II is only seen 
relatively late, and thus has a rather flat evolution with time. 
Hopefully our work inspires more detailed modeling of SN 
2009ig in the future to test our conclusions. 

These comparisons show how important it is to have the 
earliest observations possible. Out the events we consider, 
SN 201 lfe is the best constrained because it shows the largest 
velocity gradients. Just a few velocity measurements very 
early in the lightcurve can be more helpful in determining 
the explosion time than having many measurements at later 
times. Furthermore, since AMdjg oc t L76 , having observations 
only a day or two earlier probe much shallower depths in the 
ejecta. Although not discussed much here, having one or two 
early spectra that can be used for modeling the surface tem- 
perat ure can also provide tight constraints on the explosion 
time dPiro & Nak ar 2012). 

With just these three events, we are already beginning to 
see correlations between the various features that determine 
the early lightcurve rise. In the future, studies should look for 
connections between the early rise and a larger range of prop- 
erties, such as the late nebular features or the characteristics 
of the host galaxies. It will also be useful to compare spectral 
modeling methods for measuring surface abundances (like in 
lHachinger et al.l2012l) with the techniques we present here. If 
used together, they may be more constraining on the nature 
of the progenitors and the details of the explosive burning. 
Finally, it would worth exploring the early light curves of non- 
standard SNe la, like SN 2002cx dLi et al.l2003l) . Such studies 
will be important for fully utilizing the observations available 
in this new era of early detections of exploding WDs. 
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